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Abstract. The semi- classical theory of collisional depolarization of spectral lines 
of neutral atoms by atomic hydrogen (Derouich et al. 2003a Derouich et al. 2003b 
Derouich et al. 120041 and references therein) is extended to spectral lines of singly 
ionised atoms. In this paper we apply our general method to the particular cases 
of the 3d 2 D and 4p 2 P states of the Call ion and to the 5p 2 P state of the 
SrII ion. Analytical expressions of all rates as a function of local temperature are 
given. Our results for the Call ion are compared to recent quantum chemistry 
calculations. A discussion of our results is presented. 
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1. Introduction 

Observations of the linearly polarized radiation at the limb of the Sun (known as the 
"second solar spectrum"), which is formed by coherent scattering processes, show rich 
structures (Stenflo & Keller ITM71 Stenflo et al. I2UUU1 Stenflo I3UUT1 Gandorfer I2UUU1 
Gandorfer 2002 ) . The linear polarization observations reported in the atlas recently pub- 
lished by Gandorfer $20001 12002|l show significant polarization peaks in many spectral 
lines of ions, e.g. Ndll 5249 A, EuII 4129 A, Cell 4062 A, Cell 4083 A, Ba II D2 4554 
A, Zrll 5350 A, etc.. Several surveys of the scattering polarization throughout the solar 
spectrum (Stenflo et al. 119801 IT§ 83a 1983^1 see also the Q/I observations of Stenflo et 
al. 2000]and the full Stokes- vector observations of Dittmann et al. 20*t]T|l have shown that 
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ionised lines such as SrII 4078 A and the IR triplet of Call are two of the more strongly 
polarized. The interpretation of these observations leads to the indirect determination of 
the turbulent magnetic field strength via the Hanle effect. Such interpretation requires 
the solution of the coupling between the polarized radiative transfer equations (RTE) 
and the statistical equilibrium equations (SEE) taking into account the contributions of 
isotropic depolarizing collisions with neutral hydrogen. Depolarization and polarization 
transfer rates are currently available for ionised calcium levels, which have been obtained 
through sophisticated quantum chemistry methods which are accurate but cumbersome. 
Indeed, it is very difficult and sometimes not accurate to treat collision processes, in- 
volving heavy ionised atoms like Ti II, Ce II, Fe II, Cr II, Ba II..., by standard quantum 
chemistry methods. It would be useful to develop alternative methods capable of giving 
results for many levels of ionised atoms rapidly and with reasonable accuracy. 

The aim of this paper is to extend the semi-classical theory of collisional depolarization 
of spectral lines of neutral atoms by atomic hydrogen given in previous papers of this 
series (Derouich et al. 2003a Derouich et al. 2003b; Derouich et al. 2004, hereafter Papers 

1, II and III respectively) to spectral lines of ions. This paper outlines the necessary 
adjustments to the theory presented in Papers I, II and III for extension to spectral 
lines of ions. Our results are presented and compared with those obtained in the case of 
Call levels by the quantum chemistry method (Kerkeni et al. 2003). We also compare 
the results to the depolarization rates computed with the Van der Waals potential. An 
advantange of the present method is that it is not specific for a given perturbed ion, and 
may be easily applied to any singly ionised species. Indeed, we have applied our method 
to calculate depolarization and polarization transfer rates for the upper state bp 2 P of 
the SrII 4078 A line. 

The main feature of the technique is the use of perturbation theory in calculating 
the interatomic potentials. A key parameter in this theory is E p which approximates the 
energy denominator in the second-order interaction terms by an average value (Paper I 
and ABO papers: Anstee HMfl Anstee & O'Mara HMTlUMol Barklem 113551 Barklem & 
O'Mara lTMTl Barklem, O'Mara & Ross[H98j. A discussion of the effect of E p variation 
on depolarization rates is presented. Finally, we show that the present semi-classical 
method gives results in agreement with accurate but time consuming quantum chemistry 
calculations to better than 15 % for the Call ion (T= 5000K). Using our method it should 
now be possible to rapidly obtain the data needed to interpret quantitatively the Stokes 
parameters of the observed lines. 

2. Statement of the problem 

Under typical conditions of formation of observed lines in the solar atmosphere, the 
atomic system (atom, ion or molecule) suffers isotropic collisions with hydrogen atoms 
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of the medium before it radiates. The states of the bath of hydrogen atoms are unper- 
turbed. In the tensorial formulation (Fano & Racah 119591 Messiah 1961] Fano ll96*5|) . the 
internal states of the perturbed particles (here these particles are singly ionised atoms) 
are described by the spherical tensor components p k of the density matrix. Owing 
to the isotropy of the depolarizing collisions, the depolarization rates, polarization and 
population transfer rates are g-independent. The term corresponding to the depolarizing 
collisions in the master equation is given by 

J nlj Jz 

{^-^Ul = ~D k (nlJ,T)^p k 

+ Dk (nlJ -> nlJ', T) nlJ ' p k (1) 

+ quenching term. 

D k (nlJ, T) is the collisional depolarization rate of the ionic level (nlJ) at the local 
temperature T (0 < k < 2J). D°(nlJ, T) is the destruction rate of population which 
is zero since elastic collisions (J — J') do not alter the population of the level (nlJ). 
D 1 (nlJ 1 T) is the destruction rate of orientation (related to circular polarization) and 
D 2 (nlJ, T) is the destruction rate of alignment of the level (nlJ) which is of interest in 
our astrophysics framework because it is related to the observed linear polarization. 

D k (nlJ — » nlJ', T) is the polarization transfer rate between the levels \nl J) — > \nlJ'), 
where < k < k m3bX , kmax = 2 J if J < J' ( or if J > J' then fc max = 2J'). In 
particular, D k (nlJ — > nlJ', T) corresponds to collisional transfer of population (k = 0), 
orientation (k — 1) and alignment (k = 2). 

Higher order terms of D k (nlJ, T) and D k {nlJ — ► nlJ', T) with k > 2 can play 
a role in the SEE and have to be calculated. Note that, for the analysis of the linear 
polarization spectrum, only depolarization and polarization transfer rates with even k 
are need. Odd fc-terms can be eliminated from the SEE. 

D k {nlJ, T) and D k {nlJ -> nlJ', T) can be written as a linear combination of the 
collisional transition rates between the fine stucture sublevels (^{nlJMj — > nlJ'M'j, T) 
(Papers I, II and III, Sahal-Brechot I1977J) : for depolarization rates D k (nlJ, T) and 
transfer rate of population D°(nlJ — * nlJ', T), the coefficients of this linear combination 
are positive while the signs of the coefficients of the linear combination for transfer rates 
of rank k > 1 may be either positive or negative. This explains why transfer rates of rank 
k > 1 arc significantly smaller (Paper III). In our semi-classical theory, the collisional 
transition rate between the sublevels \nlJMj) — > \nlJMj) is given by (Paper I; Paper 
II): 

((nlJMj -> nlJM'j, T) = n H / 2nbdbv f(v) dv 

Jo Jo 

x\(nlJMj\I - S{b,v)\nlJM'j)\ 2 , (2) 

where f(v) is the Maxwell distribution of velocities for the local temperature T and n# 
is the local hydrogen atom number density. / is the unit matrix and T = I — S is the 
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so-called transition matrix depending on the impact-parameter b and relative velocity 
v. The collisional depolarization rates and the collisional transfer rates, which are linear 
combinations of the ((nlJMj — > nlJM'j, T) given by equation can be expressed in 
terms of the T-matrix elements. The transition matrix T is functionally dependent on 
the interaction energy matrix of hydrogen in its ground state with the perturbed ion. 
Indeed, the transition matrix elements in the dyadic basis are obtained by solving the 
time-dependent Schrodinger equation (Paper I) 

dWt)) 

(H A + H P + V eS (R))\m)=i J ^- (3) 

V e Q is the ion-hydrogen interaction used in this work and \tp(t)) is the wave function of 
the system (ion+hydrogen). Ha + Hp is the Hamiltonian of the system at the interatomic 
distance R — oo (Fig.Q. 

3. Ion-hydrogen interaction potentials 

The interaction potential for a singly ionised atom interacting with a hydrogen is treated 
in much the same way as for the neutral atom interaction with hydrogen (Papers I, II 
and III; ABO papers). In the coordinate system of Fig.^ V is given, in atomic units, by 
(Barklem & 0'Mara[ 



R + r 12 r 2 Pl ~R + r 12 r 2 Pl + ^ ind ' ' 

and atomic units are used hereafter. V- n( j = 1/R — 1/pi is the part representing an 
inductive interaction between the excess charge on the ionised atom and hydrogen atom. 
Quenching is neglected and thus we consider only the subspace nl (2/+1 states) and we 
denote the product state of the two separated atoms at R = oo by \M{). By application of 
time-independent perturbation theory to second-order, the interaction potential matrix 
elements are given by: 

(Mi\V eS \Mi) = (Mi\V\Mi) 

M^ Ml EM ' - Em 'i 

+ y mv-v^MDiMllv-v^M^ (5) 

^ — ' Em, — Em' 

Emi are the unperturbed energy eigenvalues of the isolated atoms. The expression for the 
second-order interaction can be greatly simplified if we replace the energy denominator 
Emi — Em[i of each sum, by a fixed average energy E p and assume that for important 
separations E p (R) = E p (oo). This is the Unsold approximation (Unsold 119271 Unsold 
I1955fl . E p = —4/9 atomic units is the appropriate Unsold energy value of the part of 
interaction, V^j , between excess charge on the ionized atom and hydrogen because this 
part is exactly the same as the H-H + interaction. Indeed, Unsold l|1927f) and Dalgarno & 
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Fig. 1. The perturbed ion core (with charge Z—2) is located at A and the hydrogen per- 
turbing core (a proton) at P. Their valence electrons are denoted by 1 and 2 respectively. 

Lewis (Dalgarn o fc Lewis 1956| equation 16) showed that E p = —4/9 for the long-range 
H-H + interaction. For the part of the interaction describing the interaction between the 
ion without the excess charge and hydrogen atom, the Unsold value of —4/9 cannot be 
expected to be a good approximation (Barklem & O 'Mara 1998). The reason that a value 
of —4/9 works well for neutrals is the fact that the separations of energy levels of the 
perturbed neutral atom are small compared to the separations between the ground level 
and the excited levels of the hydrogen atom, and thus the denominators are dominated 
by contributions from the H energy levels. For ions this is not the case. As a result of 
the increased core charge, the energy level spacings are generally much larger than for 
neutrals. It necessary therefore to determine E p directly for each state of the ion. The 
appropriate value of E p can be found via: 

E P = — — , (6) 

where Cq is the Van der Waals constant Cg averaged over all m substates. The Cq 
coefficient is given by the standard expression (see for example, Goodisman 1973): 

(~i 3 ^ fkk' fw ( 7 \ 

*j&fe&2+E*-E*-E*){E*-E*){E*-E*) W 

fifi and f&, are the dipole oscillator strengths of all transitions to the state of interest I for 
the perturbed ion and the ground state k for the neutral hydrogen atom. E H and E A are 
the energy eigenvalues of the hydrogen and ionised atom respectively. More details about 
the calculation of Ce are given in Barklem & O'Mara l|1998J) and references therein. The 
quantity is the mean square distance between the valence electron and the perturbed 
ion core located at A (Fig. , 

(P2>=/ P^(P2) Pi dP2, (8) 
JO 

P n 'i are the the radial wavefunctions (note P n *i{p2) — Rn*i(j>2)P2) of the valence electron 
of the perturbed atom (Anstee [19921 Seaton[1958). n* is the effective principal quantum 
number corresponding to the state \nl) of the valence electron (Papers I, II, III). 
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Using the Unsold approximation the expression for V e g becomes 

(M ; |U efr |M ; ) = (M l \V\M l )-±((M l \V\M l )) 2 + -± r {M l \V 2 \Mi) 

+ ^m\Vind\ M i)) 2 - ^(MiKdlMi)- (9) 

- j[< M ^ndl M <> ~ ((Ml^indl^)) 2 ] 

V e g of equation @ is the so-called Rayleigh-Schrodinger-Unsold (RSU) potential (ABO 
Papers). For computing V q q it is essential to determine E p in an independent calculation, 
as seen in Barklem & O'Mara l|1998l 12000). Thus for ionized atoms it is not possible 
to tabulate cross-sections as for neutral atoms (Papers I, II, III). Any calculations for 
depolarization and transfer of polarization involving ions must proceed line by line. 

4. Determination of depolarization and polarization transfer rates 

Considering a collision between a perturbed ion A and hydrogen atom H (Fig. 
Calculation of the depolarization and transfer rates follows essentially the steps listed 
below: 

1. calculation of the required atomic wavefunctions of the system A+H 

2. determination of E p directly for each state of the ion using equation © 

3. numerical evaluation of the RSU interaction energy of the system A+H given by 
equation © 

4. use of these interaction potentials in the Schrodinger equation describing the evolution 
of A+H collsional system in order to obtain the probabilities of depolarization and 
polarization transfer for a given impact parameter and a relative velocity (more details 
in Paper I; see also Papers II, III) 

5. calculation of depolarization and polarization transfer cross-sections for each relative 
velocity by integration over impact parameters 

6. integration of cross-sections over the Maxwell distribution of velocities to obtain the 
semi-classical depolarization and polarization transfer rates for a range of local tem- 
peratures of the medium 

5. Depolarization and polarization transfer rates for Ca + -H system 

An important point to emphasise is that this semi-classical method for the calculation of 
depolarization and polarization transfer rates is not specific for a given perturbed atom 
or ion. This method can be applied for any perturbed ion, but we must calculate the E p 
value for each case (section 3). Let us consider the case of the Ca + -H system in view of its 
importance in astrophysics and because it is possible to compare with recent calculations 
employing the quantum chemistry approach (Kerkenil2Q03j). The term levels associated 
to the the IR triplet lines of Call (8498 A, 8542 A, and 8662 A) are 4p 2 P 1/2 , 4p 2 P 3/2 , 
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Fig. 2. Partial Grotrian diagram of Call showing the levels and the spectral wavelengths 
in A of interest in this study. Note that the level spacings in not to scale. 



State 


(pi) (a.u.) 


C& (a.u.) 


E p (a.u.) 


3d 


7.54 


12.2 


-1.236 


4p 


22.25 


81.8 


-0.544 



Table 1. Average energy E p for the interaction of Call 3d and Ap states with hydrogen 
in its ground state together with (p\) and Cq values. 

3d 2 D 3/2 and 3d 2 D 5/2 (Figure 01. The H and K lines occur at 3969 A and 3933 A 
(Figure |2J); their upper states are also the upper states of the IR triplet. Table 2] lists, 
for the states of interest in this work, (p 2 ), C§ and the corresponding E p calculated via 
equation © (see Barklem & CMara HMSjl . 

5.1. Depolarization rates 

The depolarization transition probability is given by (Paper I; Sahal-Brechot lT^77| : 

(P k (nlJ, b, v)U = 2JTT ^ 1 {U 1 J M|7> 1 J ^ |2 
- (nlJ n\T\nlJ ii'){nl J v\T\nlJ v'Y (10) 
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' f J J k\ ( J J k \ 

X ^(_l)2./+fc+M-M 

x \ ~ v ' M' X J \ v ~M -X J 

Owing to the selection rules for the 3j-coefficients, the summation over x is reduced to 
a single term, since \ = — [p! — z/) = —(fx — v). Integration over the impact-parameter 
b and the velocity distribution for a temperature T of the medium can be performed to 
obtain the depolarization rate which is given by: 

D k {nlJ, T) ~n H vf(v)dv ^Trfog + 2tt (P k (n U,b,v)) av b dbj (11) 

where bo is the cutoff impact-parameter and we use bo = 3ao as in Anstee & O'Mara 
(1991). The excited state 4p 2 P\/2 corresponds to total angular momentum J = 1/2, the 
only non-zero depolarization rate is D 1 (4 1 1/2, T). Figure shows D 1 (A 1 1/2, T) as a 
function of the local temperature T. The non-zero depolarization rates for the 4p 2 P^/2 




L , I I I I I I I I I I I I I I I I I I I , 

2000 4000 6000 8000 1 10 4 1.2 10 4 

T(K) 

Fig. 3. Destruction rate of orientation per unit H-atom density for the Call ion, 
D 1 (4 1 1/2, T)/n H , as a function of temperature T. Z? 1 (4 1 1/2, T)/n H is given in 
rad. m 3 s . 

state are D 1 (4 1 3/2, T), D 2 (4 1 3/2, T) and D 3 (4 1 3/2, T), and these rates are 
displayed in Figure^] The non-zero depolarization rates associated to the 3d 2 D 3 / 2 and 
3d 2 D 5/2 states are D 1 ^ 2 3/2, T), D 2 (3 2 3/2, T), L> 3 (3 2 3/2, T) for 3d 2 D 3/2 and 
D 1 (3 2 5/2, T), D 2 {3 2 5/2, T), L> 3 (3 2 5/2, T), L> 4 (3 2 5/2, T) and L> 5 (3 2 5/2, T) for 
3d 2 D 5 / 2 (see Figures and 0. 



All of the rates for the 4p 2 Pi/ 2 , 4p 2 i-3/2, 3d 2 £*3/2 and 3d 2 D 5 / 2 states of Call are 
found to increase with temperature in the range 100 < T < 10000 K. As for neutral atoms, 
a functional form D(T) = BT^ X ^ 2 can usually be accurately fit to these depolarization 
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Fig. 4. Depolarization rates per unit H-atom density for the Call ion, D fe (4 1 3/2, T)/uh 
(fc=l, 2, and 3), as a function of temperature T. D k {A 1 3/2, T)/uh are given in 
rad. m 3 s . 
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1.5 10 



1 10 



5 10 



Ca + rD ) 

V 3/2' 



D 2 (2 3/2) 



(2 3/ 2) 




D 3 (2 3/2) 



2000 4000 6000 8000 1 10 4 1.2 10" 

T(K) 



Fig. 5. Depolarization rates per unit H-atom density for the Call ion, D k (3 2 3/2, T)/uh 
(fc=l, 2, and 3), as a function of temperature T. D k (3 2 3/2, T)/uh are given in 
rad. m 3 s . 

rates, where A is the so-called velocity exponent (Papers I, II, III). We find the following 
analytical expressions for the depolarization rates for 100 < T < 10000 K (except for 
L> 3 (3 2 5/2, T) and L> 5 (3 2 5/2, T) which are given for 2500 <T< 10000 K): 

- CaII(4p 2 P 1/2 )-H(Ls): 

D 1 (4 1 1/2, T) = 2.4767 x 10" 15 n H{-^—) 0A33 (rad.mV 1 ). (12) 

- CaII(4p 2 F 3/2 )-H(l, 5 ): 

D x (4 1 3/2, T) = 2.7993 x 10~ 15 n H (t^) 418 (rad.mV 1 ) 



10 
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D 5 (2 5/2): 



2000 4000 6000 8000 1 10" 1.2 10" 

T(K) 

Fig. 6. Depolarization rates per unit H-atom density for the Call ion, D k (3 2 5/2, T)/nn 
(fc=l, 2, 3, 4, and 5), as a function of temperature T. D k (3 2 5/2, T)/uh are given in 
rad. m 3 s . 



D 2 (4 1 3/2, T) = 5.2034 x 1CT 15 n H (-^)° A ° 5 (rad.mV 1 ) 



T 



D%± 1 3/2, T) = 4.8807 x lO" 10 n H (t?^) (rad.m^- 1 ) 

ouuu 



CaII(3d 2 J D 3/2 )-H(ls): 



T 



D 1 (3 2 3/2, T) = 1.9904 x 10~ 15 n H (^^)°' 373 (rad.m 3 s _1 ) 
L> 2 (3 2 3/2, T) = 1.9943 x 10~ 15 n H (-^-) - 398 (rad.mV 1 ) 

OUUU 

£> 3 (3 2 3/2, T) = 1.0772 x 10~ 15 n H (^_) - 501 (rad.mV 1 ), 



CaII(3d 2 L> 5/2 )-H(L< 



5000' 



T 



D 1 (3 2 5/2, T) = 2.0535 x 10~ 15 n H (^^)° 365 (rad.mV 1 ) 
D 2 (3 2 5/2, T) = 2.1120 x 10~ 15 n H (-^-) 384 (rad.m 3 s- 1 ) 

OUUU 

L> 3 (3 2 5/2, T) = 2.3170 x 10~ 15 n H (^Fj7j)°' 384 (rad.mV 1 ) 
L> 4 (3 2 5/2, T) = 2.0127 x 10~ 15 n H (^) 407 (rad.mV 1 ) 
L> 5 (3 2 5/2, T) = 1.2187 x 10~ 15 n H (^) 0486 (rad.m^- 1 ). 



(13) 



(14) 



(15) 



5.2. Polarization transfer rates 

The collisional transfer transition probability is given by (Paper II; Sahal-Brechot [HFfTjl : 
(P k (nlJ -> nlJ',b,v))av = (nlJ ii\T\nlJ' n')(nlJ v\T\nl J' v')* (16) 

J J k \ I J' J' k 



v' -H 1 X 
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As in equation (|llf> . the polarization transfer rates D k (nlJ — > nlJ' , T) follow from 
integration over the impact parameters and the velocities with a Maxwellian distribution. 

Inelastic collisions with neutral hydrogen which leave the radiating atom in a final 
state n'V different from the initial one nl are neglected. Only the polarization transfer 
rates inside the subspace nl are taken into account. Our transfer rates between the 
levels 4p 2 P 1/2 -» 4p 2 P 3/2 , 3d 2 D 3/2 -» 3d 2 D b/2 (D k (A 1 1/2 -> 4 1 3/2, T) and 
D k l{A 1 1/2 4 1 3/2, T)) are presented in figures □ and El respectively. D 3 (3 2 3/2 -> 
3 2 5/2, T) did not obey a power law of the form B T^ x ^ 2 . However, we can provide 
the analytical expressions for the other non-zero transfer rates: 

T 



D°(4 1 1/2 
D 1 (4 1 1/2 
D°(3 2 3/2 
D X {Z 2 3/2 
L> 2 (3 2 3/2 



4 1 3/2, T) = 4.0307 x 10~ lf W(— — )°- 407 (rad.m 3 s~ 1 ) 

OUUU 

4 1 3/2, T) = -1.1464 x 10~ i5 n ff (^— )°- 3U (rad.m 3 s _1 



-1.1464 x W~ 15 n H { 
3 2 5/2, T) = 1.8061 x 10^ 15 n ff ( 



5000' 

T - 0.392 /•_„.! I 



5000 



) u ^(rad.m d s- i ) 



T 



3 2 5/2, T) = 1J177X 10- ib n H (— — ) i4Ui (rad.m 3 s- i ) 

OUUU 

3 2 5/2, T) = 8.6286 x 10- 16 n ff (-^-)°- 490 (rad.m 3 s- 1 ) 



(17) 
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Fig. 7. Population and orientation transfer rates (k=0 and k=l), per unit H-atom den- 
sity, as a function of temperature T. The rates are given in rad. m 3 s _1 . 



6. Comparisons 

It is important to notice that the depolarization rates D k (n I J, T) as defined in this work 
(eauation llljl and the relaxation rates g k (J) as defined by Kerkeni et al. (2003) (equation 
2, section 3.2) are not equivalent. Kerkeni et al. (j2D£8j defines the depolarization cross- 
section (or relaxation rate g k {J)) as the sum of two terms: the term responsible exclusively 
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Fig. 8. Polarization transfer rates per unit H-atom density, D k (3 2 3/2 
3 2 5/2, T)/tih, as a function of temperature T. The rates are given in rad. m s . 



for the depolarization of the level (nlJ) and the term corresponding to the fine structure 
transfer between the levels (nlJ) — > (nlJ'). In our definition, D k (n I J, T) is only the 
depolarization of the level (nlJ), the fine structure transfer is not included. We calculate 
separately the rates associated to fine structure transfer which are ^-independent and 
proportional to the population transfer rates D°(n I J — > n I J' , T) (equation 4 of Paper 
II). 

For example, in accordance with our definition D°(n I J, T) = but the relaxation 
rates g°(J) from Kerkeni et al. H2003[) are not necessarily zero. In order to compare to 
the Kerkeni et al. (2003) results it is essential to substract the part of g k {J) associated 
to fine structure transfer (g(J — > J') in Kerkeni et al. I2003|) . This it is nothing more 
than a difference in definitions. Nevertheless, this difference should be taken into ac- 
count when writing the SEE. We now compare our alignment depolarization rates with 
the quantum chemistry depolarization rates (Kerkeni et al. 2003) and the alignment 

depolarization rates obtained by replacing the RSU potential with the Van der Waals 
C 

potential, V = where C§ is taken from Table 1. Note that the latter rates dif- 

fer from the usual Van der Waals formula for the rates, in that while they employ the 
van der Waals potential the collision dynamics are treated using our theory, and the 
Cq value is accurately determined (whereas typically the approximation Cq — aaip^) is 
employed where an is the polarisability of hydrogen). In Figure|5|we show our alignment 
depolarization rates with quantum chemistry depolarization rates (Kerkeni et al. I2()()3|l 
and the improved Van der Waals rates. We display only the k = 2 case which is related 
to the linear depolarization (alignment). Reference to Figure |^1 shows that the Van der 
Waals potential significantly underestimates the depolarization cross section. Our results 
show rather good agreement with quantum chemistry calculations. In concrete terms, 
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the percentage errors at T=5000 K with respect to quantum chemistry depolarization 
rates are 3.9 % for D 2 {A 1 3/2, T)/n H ] 4.6 % for L> 2 (3 2 3/2, T)/n H and 7.3 % for 
D 2 (3 2 5/2, T)/nn- The errors for the other depolarization and transfer rates are simi- 
lar to the errors for the destruction rates of alignment (in general, less than 15 %). This 
similarity is expected since all rates orginate from the same the collisional processes. 

7. Discussion 

7.1. Dependence of depolarization cross-sections on interatomic separations 

To assess the sensitivity of the calculations to the accuracy of the potentials at various 
separations we make calculations with potentials where we have introduced local pertur- 
bations. The interaction potential (Is, Apo) is multiplied by a Gaussian magnification 
factor of the form (Anstee & O'Mara lTMT)) : 

G{R) = l + cxp(-(i?-i? ) 2 ). (18) 

i?o is the position where the variation of V e Q reaches its maximum value (the interaction 
potential (Is, Apo) doubles). Figure ITUl shows the depolarization cross-section calculated 
with varying Rq. It is clear that the values of Rq inducing depolarization cross-section 
enhancement confirm the fundamental result found already for neutral atoms: the inter- 
actions of importance for the depolarization cross-section (or depolarization rate) cal- 
culations are the intermediate-range interactions (Papers I, II, and III). The principal 
differences between the RSU potentials and those from quantum chemistry occur at small 
interatomic separations. It is for this reason that we obtain rather good agreement with 
quantum chemistry calculations. The Van der Waals interaction potentials are inaccu- 
rate at the intermediate region, and this explains why these potentials underestimate the 
depolarization cross sections. 

7.2. Dependence of depolarization rates on E p 

The depolarization and transfer rates for the Ap and 3d states are calculated for 
E p = —0.544 and —1.236 respectively. As a check on the sensivity of our results to 
the precision of the calculations of E p , we have calculated the destruction rate of align- 
ment D 2 (A 1 3/2, 5000K) by varying E p in equation 10. Note that when E p decreases, 
the interaction potential decreases and so D 2 (A 1 3/2, 5000K) also decreases (Figure 
111(1 . The depolarization rate shows only an extrcmly weak variation with E p . An \E p \ 
variation of 25 %, with respect to the value E p — —0.544, corresponds to a change of less 
than 5 % in the calculated depolarization rates. It should not concluded that this is a 
general property of the depolarization rates for all states of all ionised atoms. Reference 
to the Figure UTI shows a rather strong dependence of the depolarization rates on E p for 
\E p \ < 0.25. We expect, however, that the value of \E p \ is usually greater than 0.25 and 
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the depolarization rates will not be greatly affected by possible error in the value of E p 
(see Barklem & O 'Mara 115551 l3UUU|l . 

8. Application for Sr II 4078 A line 

The Sr II 4078 A line was examined by Bianda et al. (1998), who wrote "...TTie rather 
large uncertainty in the depolarizing collision rate introduces a corresponding uncertainty 
in the field- strength scale. . ." . These authors have used the traditional Van der Waals 
approach to calculate collisional rates. The 4078 A line is the resonance line of SrII: 5s 
2 S — > 5p 2 P. We have computed the depolarization and polarization transfer rates of the 
levels 5p 2 P\/2 and hp 2 P^/2 01 SrII ion. The value of E p = —0.564 for the level 5p 2 P of 
SrII was adopted from Barklem & O'Mara (2000). We applied our method to obtaining 
the rates which are shown in Figures IT2"1 and 1T31 They were found to again obey a power 
law AT( 1 " A '/ 2 and are given for 100 < T < 10000 K by: 



- SrII(5p 2 P 1/2 )-H(ls): 



D 1 {5 1 1/2, T) 


= 2.7196 x 10" 


-15 


n H 


V 5000 ; 


(rad.m 3 s 




(19) 


SrII(5p 2 P 3/2 )- 


■H(Ls): 














D 1 (5 1 3/2, T) 


= 3.1560 x 10" 


-15 


n H 


/ \0.418 

l 5000 ; 


(rad.m 3 s~ 






L> 2 (5 1 3/2, T) 


= 5.9776 x 10" 


-15 


n H 


/ \0.406 

l 5000 J 


(rad.m 3 s~ 




(20) 


L> 3 (5 1 3/2, T) 


= 5.5413 x 10" 


-15 


n H 


/ \0.410 

l 5000 J 


(rad.m 3 s~ 







Between the term levels 5p 2 Pi/2 and 5p 2 Ps/2 there are only two non-zero polarization 
transfer rates which are given in Figure IT^l The analytical expressions for these rates for 
100 < T < 10000 K are: 



D°(5 11/2^51 3/2, T) = 4.6184 x 10" 15 n H ( ) 409 (rad^s -1 ) (21) 

5000 

D 1 (5 1 1/2 -»■ 5 1 3/2, T) = -1.5309 x 10~ 15 n H (-L-f-™ ( ra d. m 3 s- 1 ) 



9. Conclusion 

We have adapted our semiclassical method of calculation of collisional depolarization of 
spectral lines of neutral atoms by atomic hydrogen to allow it to be used for singly ionised 
atoms. Comparison with recent quantum chemistry calculations for Call indicates and 
error at T=5000 K of less than 5 %. This is an encouraging result which supports the 
validity of our semi-classical approach. Using this method we should be able to calculate 
depolarization rates of the levels involved in transitions of heavy ionised atoms like SrII, 
Ti II, Ce II, Fe II, Cr II, Ball... Calculations must proceed line by line because a suitable 
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E p value needs to be determined for each relevant state of the given ion. We have applied 
the method to calculate depolarization and polarization transfer rates for the SrII 5p 2 P 
state. These calculations should allow a more accurate theoretical interpretation of the 
observed linear polarization of the SrII 4078 A line. 
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Fig. 9. Depolarization rates for k = 2 as a function of temperature. Full lines: our results; 
dotted lines: quantum chemistry calculations fKerkeni I2(J(J3[I : dot-dashed lines Van der 
Waals approximation. 
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Fig. 10. Depolarization cross section enhancement due to a Gaussian local perturba- 
tion of the potential. Cross sections are calculated for the relative velocities: 3 km s _1 , 
10 km s^ 1 and 17 km s _1 . 
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Fig. 11. Dependence on E p of the destruction rate of alignment D 2 (4 1 3/2, 5000K)/n#. 
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Fig. 12. Destruction rate of orientation per unit H-atom density for the SrII 
ion, D l (f> 1 1/2, T)/riH, as a function of the temperature of the medium, T. 
£> 1 (5 1 1/2, T)/riH is given in rad. m 3 s _1 . 
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Fig. 13. Depolarization rates, D fc (5 1 3/2, T)/nn (k=l, 2, and 3), as a function of 
temperature T. D k (5 1 3/2, T)/uh are given in rad. m 3 s _1 . 
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Fig. 14. Population and orientation transfer rates (fc=0 and k=l) for the SrII ion, per 
unit H-atom density, as a function of temperature T. The rates are given in rad. m 3 s _1 . 



